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ABSTRACT 

Smoothed particle hydrodynamics (SPH) is used to estimate accretion rates of mass, 
linear and angular momentum in a binary system where one component undergoes 
mass loss through a wind. Physical parameters are chosen such as to model the alleged 
binary precursors of barium stars, whose chemical peculiarities are believed to result 
from the accretion of the wind from a companion formerly on the asymptotic giant 
branch (AGB). The binary system modelled consists of a 3M Q AGB star (losing mass 
at a rate 10~ 6 M Q y" 1 ) and a 1.5M0 star on the main sequence, in a 3 AU circular 
orbit. Three-dimensional simulations are performed for gases with polytropic indices 
7 = 1, 1.1 and 1.5, to bracket more realistic situations that would include radiative 
cooling. Mass accretion rates are found to depend on resolution and we estimate 
typical values of 1-2% for the 7 = 1.5 case and 8% for the other models. The highest 
resolution obtained (with 400k particles) corresponds to an accretor of linear size 
ks 16Rq. Despite being (in the 7 = 1.5 case) about ten times smaller than theoretical 
estimates based on the Bondi-Hoyle prescription, the SPH accretion rates remain large 
enough to explain the pollution of barium stars. Uncertainties in the current SPH rates 
remain however, due to the simplified treatment of the wind acceleration mechanism, 
as well as to the absence of any cooling prescription and to the limited numerical 
resolution. 

Angular momentum transfer leads to significant spin up of the accretor and can 
account for the rapid rotation of HD 165141, a barium star with a young white dwarf 
companion and a rotation rate unusually large among K giants. 

In the circular orbit modelled in this paper, hydrodynamic thrust and gravita- 
tional drag almost exactly compensate and so the net transfer of linear momentum is 
nearly zero. For small but finite eccentricities and the chosen set of parameters, the 
eccentricity tends to decrease. 

Key words: accretion, accretion discs - hydrodynamics - binaries : close - stars: 
barium 



1 INTRODUCTION 

Wind accretion in a binary system influences the surface 
composition and spin of the accreting star and can change 
the orbital parameters of the binary. This paper aims at es- 
timating these effects in a binary system in which one of 
the components is an asymptotic giant branch (AGB) star 
which suffers from strong mass loss (M ~ 10 -6 Mq y _1 ) 
through a wind. The companion main sequence (MS) star 
accretes some of this carbon and s-process element enriched 
gas and mixes it in its envelope. Stellar evolution will even- 
tually transform the AGB star into a white dwarf (WD) and 
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the initially less massive MS star into a chemically peculiar 
giant. 

This is the evolutionary scenario proposed by Boffin 
& Jorissen (1988) for the formation of barium stars, i.e. 
late-type giants exhibiting overabundances in carbon and s- 
process elements (e.g. Lambert 1985). All barium stars are 
likely members of a binary system with a WD companion 
(McClure & Woodsworth 1990, Jorissen & Boffin 1992). 

The theoretical description of the accretion process is 
complex due to the intricate structure of the flow and the 
complex non-linear hydrodynamics of the shocked wind ma- 
terial. One can estimate the expected fraction M acc / M of 
material being accreted (Boffin & Jorissen 1988) from an in- 
terpolation formula due to Bondi & Hoyle (1944) and Bondi 
(1952) for the accretion rate onto a star moving at con- 
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stant velocity through a gas of uniform density and tem- 
perature. This involves the interpolation of the estimated 
accretion rates from two even more simplified models. In 
the first model (Hoyle & Lyttleton 1939), one neglects gas 
pressure and estimates the accretion rate Mhl onto a com- 
pact gravitating object with mass M, moving at constant 
velocity Voo through a pressureless medium with initially 
uniform density poo . In this model, gravitationally deflected 
material which passed on one side of the star collides with 
material passing on the opposite side, thereby cancelling its 
transverse velocity and forming an 'accretion line' behind 
the star. The material from this accretion line that has a 
velocity below the local escape velocity from the star will be 
accreted and one finds (Hoyle & Lyttleton 1939): 

Mhl = ttRaPooVoo, (1) 

where the accretion radius Ra is given by 

R A =2GM/v 2 00 . (2) 

The second model (Bondi 1952; see also Theuns & David 
1992 for the closed form solution of the flow pattern) in- 
cludes gas pressure but neglects instead the motion of the 
star through the medium. The accretion rate in a steady 
state is not uniquely determined but the model does pro- 
vide a maximum accretion rate Mb'- 

M b = ctnRspooCoo, (3) 

where the Bondi radius is given by Eq. (^) but with the 
sound speed at infinity Coo replacing Voo- The efficiency pa- 
rameter q is of the order of unity and depends on the value 7 
of the polytropic index of the gas. The accretion rate Mbh 
for the Bondi-Hoyle case, where both the velocity of the star 
as well as gas pressure are included, is based on an interpo- 
lation between these two extreme cases: 

1 + ^ 2 J , (4) 

where Moo = Uoo/coo is the Mach number. 

The theoretically predicted accretion rate Mbh in the 
Bondi-Hoyle model has been extensively tested against nu- 
merical simulations by e.g. Hunt (1971) and more recently 
by Ruffert (1994) and Ruffert & Arnett (1994). The theo- 
retical and numerical rates agree to within 10% which must 
be considered a bit fortuitous since e.g., although transiently 
present, no accretion line (as envisaged in the Hoyle and Lyt- 
tleton picture) actually exists. The latter simulations also 
show the importance of numerical resolution on the com- 
puted accretion rates. 

In a previous paper (Theuns & Jorissen 1992, Paper I), 
we presented three-dimensional (3D) smoothed particle hy- 
drodynamics (SPH) simulations of wind accretion, taking 
into account the binary motion. Two-dimensional simula- 
tions of this problem were previously done by S0rensen et 
al. (1975) and by Matsuda et al. (1987, 1992). Even a glanc- 
ing comparison between the structure of the flow in this 
case (Paper I) against that in the Bondi-Hoyle model (as in 
e.g. Ruffert & Arnett 1994) shows that the binary motion 
induces a very different flow pattern. Consequently, using 
the accretion rate Mbh in this situation as well might be 
stretching one's luck too far. 

The main aim of this paper is to test whether the mass 
accretion rates are large enough to produce barium stars 



through the wind accretion scenario discussed above. The 
impact of wind accretion on the orbital parameters and on 
the spin velocity of the accretor will also be evaluated. An es- 
pecially important diagnostic of wind accretion is provided 
by the variation of the orbital eccentricity caused by the 
transfers of linear momentum between the gas and the ac- 
creting star. Moreover, if the accreted gas possesses angular 
momentum, the spin period of the accreting star can be al- 
tered. Several aspects considered in the present paper may 
also be relevant to the study of wind accretion in symbiotic 
stars or in C, Aurigae systems. 

In addition to the complications introduced by the bi- 
nary motion, uncertainties about the thermodynamic state 
of the gas call for caution when deducing the accretion rate 
from oversimplified models. In fact, the simple polytropic 
equation of state adopted in this paper for the gas may not 
be appropriate in these systems exhibiting 'cooling radia- 
tion' in the form of UV light. Interacting systems such as 
symbiotic or extrinsic S stars, where wind accretion is cur- 
rently taking place, indeed exhibit a wealth of UV emission 
lines or continuum UV radiation (e.g. Nussbaumer & Sten- 
cel 1987; Johnson & Ameen 1991). In the interacting S star 
HD 35155 for example, a total UV continuum emission of 
about 0.2 L© is observed (Johnson & Ameen 1991). In Pa- 
per I, we identified regions of shocked gas as being respon- 
sible for this kind of emission. 

Comparison of the accretion rates obtained for different 
7's nevertheless allows us to appraise the importance of the 
assumed polytropic index on the computed rates. 

2 MODELING WIND ACCRETION WITH SPH 

2.1 Description of the numerical method 

All numerical details can be found in Paper I. Only a short 
summary will be given here. Calculations are done in an 
inertial frame in which the motion of the centers of mass of 
the two stars, AGB and MS, are computed from the solution 
of the two-body problem for a circular orbit, neglecting all 
hydrodynamic effects. The hydrodynamics of the gas, on the 
other hand, is treated by means of variable resolution SPH 
(Lucy 1977; Gingold & Monaghan 1977; see e.g. Benz 1989 
for a review). 

Gas forces include the hydrodynamic pressure force, 
Vp/p, the gravitational force from the MS star and from 
the AGB star. In addition, a mechanism responsible for the 
AGB wind acceleration is included, simulating the effects of 
the AGB pulsations and radiation pressure on the escaping 
gas (see e.g. Bowen 1988). Since the luminosity force per 
unit mass falls with distance like 1/r 2 , just as the gravi- 
tational force does, one can implement such an acceleration 
mechanism by decreasing the effective AGB star mass where 
it enters the interaction with the gas. In these simulations, 
we assume an accelerating force which exactly balances the 
gravitational force from the AGB (resulting in a wind veloc- 
ity independent of the distance from the mass-losing star, as 
obtained in Bowen's wind models) . Self-gravity of the gas is 
neglected. 

A specific model is characterized by the orbital elements 
of the binary (radius A of the circular orbit and period P 
corresponding to a circular velocity V c = 2-nA/P), the mass 
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ratio of the two stars and the properties of the wind, i.e. 
its polytropic index 7 and its velocity, mass flux and sound 
speed at the surface of the AGB star. 

In addition to these physical parameters, there are addi- 
tional purely numerical parameters, most critically the effec- 
tive resolution. The latter is determined by the total num- 
ber of SPH particles used. We estimate this resolution by 
computing the average smoothing length h (see Paper I) for 
particles around the accreting star. That length character- 
izes the size of the interpolation region over which quantities 
are smoothed in the SPH formalism. 

The simulations are started with vacuum initial condi- 
tions. A supersonic wind is started around the AGB star and 
the model is integrated until a quasi steady-state is reached. 
The resolution is then improved by increasing the number of 
particles emanating from the AGB star per unit time (yet 
keeping the mass flux fixed). It is possible to monitor in 
this way how the accretion rates change as a function of 
resolution. The accretion process is simulated by removing 
mass from the flow, once it enters a sphere with radius h 
around the accreting star. The average resolution h is com- 
puted from all particles in a sphere of radius th around the 
accreting star. 

2.2 Calculation of fluxes 

Accretion rates are computed by measuring fluxes round the 
secondary. Two different methods are used to compute fluxes 
(see also Paper I). The flux <&a of a quantity A is defined 
by 



4> 



Apv ■ dS, 



(•») 



where p and v denote the fluid density and velocity respec- 
tively, dS is the outward surface element of a sphere with 
radius R centered on the accreting star. In the first method 
for evaluating fluxes, the flux density Apv entering Eq. Q 
has been evaluated using its SPH estimate on 20x20 points 
uniformly spaced in cos 9 and <j> on the surface of the sphere 
of radius R (with 9 and <f> being the usual spherical coordi- 
nates with respect to the centre of the sphere). A numerical 
integration over the surface of the sphere then gives &a- 

Alternatively, &a can be computed by summing over 
SPH particles of mass m^: 



,, \ -> dim 



(6) 



For A = 1, &a represents the mass accretion rate. In 
that case, <f>A should be independent of R for a stationary 
flow (as long as the sphere does not reach other sources or 
sinks, like e.g. the mass-losing star or the outer boundary). 
The comparison of the mass accretion rates computed from 
the flux through spheres of different radii thus provides a 
way of evaluating the stationarity of the flow (see Fig. |^). 
For A = v x , &a represents the momentum accreted in the 
x-direction. In that case, &a does depend on R, since the 
velocity v x of a fluid particle changes due to the gravitational 
force of the accreting star. However, the total force on the 
star, which is the sum of the accreted momentum and the 
gravitational force exerted by the gas on the star, should 



2.3 Test simulations 

We tested our code and the prescription for calculating 
fluxes by simulating Bondi accretion, i.e. the spherically 
symmetric accretion of polytropic gas onto an accretor of 
constant mass M (Sect, [fj) at rest. Given the sound speed Coo 
and density p^ of the gas at infinity, this problem has a fam- 
ily of solutions which differ by the effective accretion rate. 
The so-called transonic branch corresponds to the situation 
where the flow changes from subsonic at large distances to 
supersonic close to the accreting object. This solution has 
the largest accretion rate (e.g. Theuns & David 1992). A 
flow pattern corresponding to this branch and matching the 
analytical solution is first set up. The parameters M, Coo and 
poo are then scaled in order to get flows with the same aver- 
age resolution h within the sphere of radius Th, but with the 
sonic point falling either outside the accreting sphere's sur- 
face (the 'supersonic' model) or inside the accreting sphere 
(the 'subsonic' model). Note that the physical model is iden- 
tical in both cases, the only difference being in the place 
where mass is removed from the flow: in the subsonic model, 
potentially all mass is removed from the flow before it has 
managed to become supersonic. 

The upper panel of Fig. [l] compares the simulation re- 
sults for the 'subsonic' case with the analytical solution. The 
agreement is satisfactory outside the accreting sphere (only 
there is the comparison meaningful, obviously). The deduced 
accretion rates in the stationary regime are 20% low with re- 
spect to the analytic value (lower panel of Fig. jjj) , which is 
sufficient for our purpose. The agreement is better in case 
the sonic point falls outside the accreting sphere (the 'su- 
personic' model). Note that the difference between analytic 
solution and simulation is likely to be due, not only to the 
accuracy of the flux calculation itself, but also to how well 
the initial profile matches the analytic profile. The resolu- 
tion in these test cases, in terms of the number of particles 
inside the region in which particles are being accreted, is 
comparable to the resolution obtained in the standard wind 
accretion case, to be discussed in Sect. 2.4. 



again be independent of R. Finally, for A - 
is the accreted spin momentum. 



yv x , $a 



2.4 Numerical simulations 

All the simulations considered in this paper are for a bi- 
nary system consisting of a 3 Mq AGB star and a 1.5 Mq 
companion in a circular orbit around the AGB star with a 
semi-major axis of 3 AU (corresponding to a period of 895 d, 
and an orbital velocity of 36 km s _1 ). The AGB star is losing 
mass at a steady rate of Mi = 10 -6 Mq y _1 with a (con- 
stant) wind velocity of 15 km s _1 (here, and in what follows, 
subscripts 1 will refer to the mass-losing star). The sound 
speed at the base of the flow is 7.6 km s _1 (corresponding to 
a temperature of about 4050/7 K at the surface of the AGB 
star; note that this sound speed is larger than the value used 
in Paper I). Three models are considered which differ only 
in the assumed value of 7 (1.5, 1.1 and 1). Unless stated 
otherwise, we use units in which period P, semi-major axis 
A and gravitational co nsta nt are all unity. In these units, 
the parameter Th (Sect. 2.1) is taken equal to 0.3. 

Table |l| provides the conditions for the different runs 
performed and the corresponding accretion rates. The suc- 
cessive columns give the total number of particles N, the 
polytropic index 7, the average resolution 2h inside the 
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Figure 1. Comparison of SPH test simulation (with 7 = 1.5) against analytic solution for the transonic branch of Bondi accretion. The 
corresponding Bondi radius, as defined in Eq.fel), is 400 in our units. Upper panel: SPH profiles (in the stationary regime) of density 
p, infall velocity v r , sound speed c, mass flux A = —Airr 2 pv r (normalised by the analytical value Ag) and Mach number M compared 
with the analytic solution (solid lines). For the sake of clarity, density, velocity and sound speed profiles have been offset by -1, -0.7 and 
-0.4, respectively. Particles to the left of the vertical line are used to compute the average resolution h around the accreting object (i.e. 
l°gio(^) = 0-9)- 

Lower panel: accretion rates (normalised to the analytical value) determined from Eq. (ph (curve labelled 'acc') compared with rates 
determined from Eq. (0) for three values of the radius R (R = 2h, labelled 'g2'; R = 3h, labelled 'g3'; R = 4h, labelled 'g4', where 
h 0.8), as a function of time. For times > 2, a steady-state is reached. 



Table 1. Basic data and accretion rates for the different cases considered 



N 


7 


2h/A 


2h 
R 


t/P 


Ms /Mi 

% 


Hy 

M 2 V C 


Gy 

M 2 V C 


Fy 
M 2 V C 


L Z /V C AM 2 


100k 


1.5 


0.11 


71 


24 


2 


-1.5 


1.55 


0.05 


0.06 


100k 


1.1 


0.10 


61 


17 


8 








0.06 


70k 


1 


0.10 


64 


22 


8 








0.06 



sphere of radius th around the accreting star [as is common 
practice, we use 2h to characterise the resolution because 
the SPH kernel W(r/h) becomes zero for r > 2h, see Pa- 
per I], expressed relative to the orbital separation and in 
solar radii, the number of orbital periods during which the 
simulation has been carried out (t/P), the fraction of the 
mass lost by the AGB star that is accreted by the companion 
(M 2 /Mi), the hydrodynamic, gravitational and total drags 
(H y , Gy and F y , respectively, expressed in terms of the ac- 
creted momentum M 2 V C ), and finally the dimensionless spin 
accretion rate L Z /V C AM 2 . The various accretion rates will 
be described in more details in Sect. 3.1. A dash in Table [l] 
means that the value has not been computed. 

Each simulation typically requires about two weeks of 
CPU time on a HP-730 workstation. In terms of the accre- 
tion radius Ra [Eq.(^)], the resolution achieved around the 
accreting star amounts to 2R/Ra ~ 0.2, since Ra = 380 R©. 
Since this resolution is much higher than in Paper I (100k 
particles instead of 40k particles for the case 7 = 1.5), Fig. [| 
displays a general view of the flow pattern, as well as a zoom 
on the region close to the accreting star. Despite a larger 
sound speed at the surface of the AGB star, the general 



structure of the flow is similar to that described in Paper I 



5.1 



It is worth noting in relation with the discussion of Sect 
that the Roche-lobe geometry is not at all apparent in the 
flow structure, because the wind acceleration mechanism has 
the net effect of suppre ssin g the gravitational potential of 
the AGB star (see Sect. Et|). 



3 RESULTS 

3.1 SPH accretion rates 

The accretion rates of mass, momentum and spin for the 
three values of the polytropic index (7 = 1,1.1 and 1.5) 
considered in this paper are given in Table [l]. They will be 
discussed in some detail here for the 7 = 1.5 (adiabatic) 
case only, since the same remarks apply to the other cases 
as well. 

Figure |^ shows the accretion rates as a function of time. 
The top panel indicates that the two methods of determin- 
ing the fluxes, namely direct summation over SPH particles 
(Eq. ^ or estimates of SPH fluxes (Eq. ^) , give very similar 
results: about 2 to 2.5% of the mass lost by the AGB star is 
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Figure 2. Left panel: General structure of the velocity field in the orbital plane for the adiabatic case 7 = 1.5 with 100k particles. Only 
particles lying in a z-slice [-0.1,0.1] enclosing the orbital plane are represented. The stationary frame in which the mass-losing star is at 
(x = 0.33, y = 0, z = 0) and the accreting star at (x = —0.66, y = 0, z = 0) has been used. 

Right panel: Same as left panel, zoomed in on the accreting star. The spin momentum carried by the wind close to the accretor is 
apparent. 



accreted by the companion (but see below) . Incidentally, the 
fact that the mass fluxes through spheres of different radii 
are almost equal indicates that the flow pattern is relatively 
stationary. 

The middle panel gives the forces exerted by the gas on 
the accreting star, due to the accretion of linear momentum 
(hydrodynamic drag, labeled H) and due to the gravita- 
tional pull of the gas lagging behind the accretor (labeled 
G). More precisely, the accretion of linear momentum has 
been computed from 

H = P-Af 2 v 2 , (7) 

where P e ^ v w dm w /dt denotes the momentum in the 
wind, accreted by the star. The index w indicates that the 
sum extends over SPH particles. The inertial forces in this 
figure are projected on a set of axes, with the x axis along 
the line joining the two stars, the mass-losing star lying on 
the +x axis with its orbital motion pointing along +y, and 
finally the z axis being perpendicular to the orbital plane. 
The gravitational drag G, on the other hand, is computed 
by summing over SPH particles in the wind as: 

G = ^ Fg r av,w- (8) 

w 

The net force acting on the star is labeled F = G + H . 

The hydrodynamic drag is negative, for both radial and 
tangential directions. This means that the accreting star is 
propelled away from its companion and is accelerated along 
its orbit due to the gas raining back onto it from y > 
(i.e., a thrust rather than a drag), as apparent from Fig. || 
However, the gravitational forces have exactly the opposite 



tendency and the net effect is an inward pull and a very 
slight slowing down of the star along its orbit. 

The bottom panel shows the accretion of the dimension- 
less spin angular momentum, again obtained by summing 
over SPH particles: 

— => (r„ - r 2 ) x (v„ - v 2 ) i — (9) 

AV C M 2 ^ AV C M 2 dt 

w 

The radial (a;) and tangential (y) components fluctuate 
around zero, as expected from the symmetry of the flow 
along the z direction (perpendicular to the orbital plane). 
There is, however, a net accretion of spin momentum in the 
z direction, in such a way that the accreting star tends to be 
spun up if it is in synchronous rotation to begin with. This 
can be clearly seen on Fig. |^ from the direction of rotation of 
the disc-like structure around the accreting star. The accre- 
tion rates of dimensionless spin momentum are remarkably 
similar amongst the three models, notwithstanding the fact 
that a much better-developed accretion disc forms in the 
isothermal case (compare Figs. 1 and 8 of Paper I). 

We tried to improve our estimate of M2/M1 by in- 
creasing the resolution drastically in a small region around 
the accreting star. However, conclusions drawn from sim- 
ulations involving such a restricted region will be mean- 
ingless once boundary effects propagating inwards become 
important. In order to illustrate that effect, a simulation 
has been performed with most of the particles removed, ex- 
cept for those close to the accreting star. As illustrated by 
Fig. [|, the accretion rate does not change when measured 
over a sufficiently short interval of time, but starts deviating 
when boundary effects propagate into the accretion sphere 
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(round t/P = 24.34). In Fig. §, the line labeled 'large box' 

corresponds to the accretion rate for the standard modelQ 
whereas the dotted line, labeled 'small box', includes parti- 
cles close to the accreting star only (in fact, the small box 
with the same resolution as the large box, contains only 4% 
of the number of particles of the large box). Obviously, for 
t/P < 24.34, the accretion rates remain identical yet start 
to diverge afterwards. Therefore, the high-resolution simula- 
tions described below, and involving the small box, are only 
meaningful for t/P < 24.34. In order to artificially increase 
the resolution, we proceed as follows. At time t/P = 24.30, 
each SPH particle is replaced by a swarm of n particles, 
distributed around their mother according to the SPH den- 
sity kernel. This increases the resolution by a factor n 1 ^ 3 . 
The different curves in Fig. ^ correspond to simulations with 
n = 1 (dotted), n = 10 (triangles), n — 25 (dots) and n — 90 
(squares) , which translates into a number of particles in the 
small box around the accreting star of 4k, 40k, lOOfc and 
360fc, respectively. Drawing conclusions from these numeri- 
cal experiments is hampered by the fact that the accretion 
rate varies quite considerably over time intervals of the or- 
der of 0.02P, as is clear from the n = 1 case (solid line). 
These variations on a time scale much longer than the time 
step St < 5 x 10 _4 P simply reflect the discrete episodes 
of particle ejections by the mass-losing star (see Sect 2.3 of 
Paper I). In addition, once SPH particles are cloned, as de- 
scribed above, the system needs some time to relax to a new 
steady state. The following conclusions can nevertheless be 
drawn from these simulations: 1) the accretion rate decreases 
with increasing resolution, and 2) the mass accretion rate is 
not likely to be much less than about 1%, i.e. about 1/3 lower 
than deduced from the standard model. This value follows 
assuming that the almost linear trend observed in Fig. ^| 
which displays numerical resolution versus mass accretion 
rate, can be extrapolated to the size of a main sequence ac- 
cretor (i.e., 2h ~ 1.5 x 10~ 3 which corresponds to IRq). In 
what follows we will use an estimated mass accretion rate 
of 1% for the adiabatic model, but we stress that an accu- 
rate estimate of the real accretion rate requires much higher 
resolution simulations than the ones presented here. 



3.2 



Comparison with the Bondi-Hoyle mass 
accretion rate 



In Sect. 3.1, a mass accretion rate /3 = —dM^/dMx of the 
order 1% is derived from our SPH simulations in the 7 = 1.5 
case. That result will now be compared to the accretion 
rate derived by the Bondi & Hoyle (1944) formula. In the 
framework of a binary system, the gas velocity v^, appearing 
in Eq. (Q) may be replaced by the relative velocity between 
the wind and the accreting star, whereas the gas density poo 
may be replaced by the wind density p at the location of the 
accreting star. The latter may be deduced from the mass 
conservation equation Mi = —4Tvr 2 v vi p. Substituting these 
expressions into Eq. (W) yields 
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Figure 5. Accretion rates of mass (M2/M1, in %) as a function 
of resolution, 2h, for the 7 = 1.5 model. The upper and lower 
curves correspond to the maximum and minimum accretion rates 
for these resolutions, read from Fig. Q 



a_ / GMi \ ' 
A 2 \ vl ) 



[l + KrbA>w) 2 + (cA; w ) 2 ] 3 / 2 ' 
or using Kepler's third law: 



A: 4 



[1+k 2 + (c/Ww) 2 ] 1 - 5 ' 



(10) 



(11) 



where a is a parameter of order unity not fixed by the theory, 
p = Mi/(M\ + M2), k = Uorb/^w, and v OTh is the orbital 
velocity. 

For the conditions given in Sect. 2.4, Eq. (|ll|) yields 
(3 ~ 10 to 20% for k = 2.4, taking a between 0.5 and 1, as 
suggested by the numerical simulations of Ruffert & Arnett 
(1994). The accretion rate deduced from the SPH simula- 
tion in the 7 = 1.5 case is thus about ten times smaller than 
that predicted by the Bondi & Hoyle formula. The reason for 
this apparent discrepancy is that Eq. (jl^) is only applicable 
to situations where k « 1, since it was originally derived 
by considering a single star moving through a gas cloud. In 
the case of a fast wind [k « 1), the concept of accretion 
column as introduced originally by Bondi & Hoyle to de- 
scribe the accretion process, does hold with the accretion 
column being tilted with respect to the radius vector. When 
k is of the order of unity, however, as is the case here, the 
binary motion strongly disturbs the shape of the accretion 
column, as is obvious from Fig. |j| The fact that the above 
formula becomes meaningless for large values of k is further 
illustrated by the fact that (5 increases like k when k >> 1, 
though values of f3 larger than unity have of course no physi- 
cal meaning. The fact that the Bondi-Hoyle formula Eq. ( pd| ) 
was not expected to hold in the case of accretion of a slow 
wind was the prime motivation to perform the simulations 
described here. 

Finally, note that the term (c/« w ) 2 in Eq. ( |ll| ) is usu- 
ally negligible (in the case of adiabatic wind expansion, c 



decreases as r 



-(7-1)1 



and has therefore not been consid- 



ered in the above discussion. 



t The reason why the standard case yields accretion rates of the 
order of 2.5-3.5% in Fig. § as compared to 2.0-2.5% in Fig. ||is 
because the latter correspond to time averages. 



4 CONSEQUENCES OF WIND ACCRETION 

Wind accretion can alter the evolution of a binary in several 
respects: 
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Figure 3. Accretion rates as a function of time for the 7 = 1.5 case. 

Top panel: accretion rates of mass (M2/M1, in %)■ The different symbols represent estimates based on Eq. ^ (filled dots, labeled 'acc') 
or on Eq. (^) with three different radii R (0.11, 0.16 and 0.22, labeled l g2', 'g3' and l g4', respectively). 

Middle panel: accretion rates of linear momentum in the radial (subscript x) and tangential (y) directions, in units of M^y c \ H refers to 
the hydrodynamic drag, G to the gravitational drag and F to the sum of those (see text for details). 

Bottom panel: rate of accretion of dimensionless spin momentum, L / M^VcA, in the radial x, tangential y and 2 directions. Spin momenta 
along x and y are multiplied by ten. 



Figure 4. Accretion rates of mass (M2/M1, in %) as a function of time (in units of the orbital period P) for the 7 = 1.5 model with 
different resolutions 2h. The solid curve labeled 'large box' refers to the standard model. The dotted line labeled 'small box' corresponds 
to a simulation where only those particles close to the accretor are retained. Up to t/P m 24.34, these two rates are identical. The other 
curves (triangles, dots and squares) refer to simulations in a 'small box' as above, but with an initial number of particles increased as 
indicated in the text. The different curves are labeled with their effective resolutions. 




1 
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• the chemical composition of the envelope of the accret- 
ing star may be modified if the accreted matter has a differ- 
ent composition, and if enough such mass is accreted; 

• the orbital elements may change as a result of the mass 
loss and mass transfer; 

• the accreting star may be spun up if the accreted matter 
possesses significant angular momentum. 

Each of these effects is now considered in turn on general 
grounds. Application to barium stars will be considered in 
Sect. B. 



4.1 Envelope pollution 

The pollution of the envelope of the accreting star depends 
on (i) the amount of mass AM 2 accreted from the wind, (ii) 
the dilution factor of the accreted matter in the envelope, 
and (iii) the extent to which the chemical composition of 
the accreted matter differs from that of the envelope. We 
will assume that the accreted matter has been fully mixed 
in the envelope. In the case that the accreting star is on the 
giant branch, such a mixing naturally results from the con- 
vective nature of the envelope. In the other case where the 
accreting star is still on the main sequence, some mixing will 
occur if the accreted matter has a mean molecular weight 
larger than that of the underlying layers (Proffitt 1989; Prof- 
fitt & Michaud 1989). Such an inversion of mean molecular 
weight, triggering a turbulent mixing until the inversion is 
suppressed, is expected if carbon-rich matter from the wind 
of an AGB star falls on top of the surface layers of the main 
sequence star. 

The overabundance /; of element i in the envelope of 
the accreting star (i.e. the ratio between the abundance af- 
ter completion of the accretion and mixing processes, and 
the abundance in the primordial envelope) is related to its 
overabundance gt in the wind (i.e. in the AGB atmosphere) 
through the relation 



giAMg + M 2 ,o - M 2 , c 
AM 2 + M 2 , - M 



2, core 



(12) 



where T is the dilution factor of the accreted matter AM 2 
in the envelope of mass (M 2 — M2, C ore)- Here, M 2t o and 
M 2tCOle denote the total initial and core mass, respectively, 
of the accreting star. In deriving the above relation, it has 
been assumed that gas and grains are accreted with the 
same efficiency. The chemical fractionation that occurs be- 
tween the gas and grain phases in the AGB wind (due to 
the different refractory properties of the different chemical 
elements) is then erased once they mix again in the com- 
panion's envelope. This assumption is actually not always 
satisfied: in some observed post- AGB star binaries only gas 
and no grains have been accreted (e.g. Waters et al. 1992). 



4.2 Variation of orbital parameters 

The orbital parameters of the binary will change as a con- 
sequence of the mass lost from the system, the mass trans- 
ferred from one star to the other, the gravitational force of 
the wind on both stars and finally, the change in momentum 
of the accretor, due to the accretion of momentum from the 
wind. In addition, the mass loss process itself may be asym- 
metric, providing an extra source of momentum changing 



the orbital parameters, but that possibility will not be con- 
sidered here. 

Huang (1956) derives expressions for such changes in an 
adiabatic approximation, i.e., when the change per period is 
small. However, Huang's parametrization does not allow to 
compare the respective effects of drag and mass accretion 
on the variations of the orbital elements, though the for- 
mer effect is of importance in the situation modelled in this 
paper (as can be judged from Fig. [|). Therefore, Eqs. (|l^ 
- |l5| ) describing the variation of the orbital elements are 
parametrized in term of F y , the tangential component (i.e. 
perpendicular to the radius-vector, and positive along the 
orbital motion of the mass-losing star 1) of the total force, 
including gravity and acc retion of linear momentum. With 
O, F v = G 2 . v - {M 2 /M l )G 1 , y 



the notations of Sect. 



+ H y 



as displayed in Fig. || (but note that we neglected the grav- 
itational force Gi, M exerted by the gas on the mass-losing 
star in that figure). That quantity, as well as the mass ac- 
cretion rate, are assumed to remain constant during one or- 
bital cycle: this is probably a good approximation as long as 
e <C 1 (e is the eccentricity of the orbit). Following this line 
of argument (see the Appendix for more details), we find, 
neglecting all terms of order e 2 or higher: 



A 
A 
P 
P 



= -2 



M 1 +M 2 
Mi + M 2 
Mi + M 2 
'M! + M 2 



M 2 V C 
'M 2 V C 



3 Fy 



1M2 

' 2M 2 2M 2 V C 



(13) 
(14) 
(15) 



Equations (|l^) and ( jl^ ) illustrate basically the same phys- 
ical effects: mass loss from the system (Mi + M 2 < 0) in- 
creases A and P (first term) , and slowing down of the accre- 
tor or speeding up of the mass-loser due to the exerted forces 
(F y > 0) has the opposite effect (last term). The behaviour 
of the eccentricity can be written slightly more revealing in 
the case of equal masses: 



3.1 
2 ( 6 + 



w,y \ M 2 
V c ' M 2 

(for Mi = M 2 ) 



3 Gi, v — G\,y 
2 M 2 V C ' 



(16) 



where v WlV is the tangential wind velocity, at the point it 
is accreted. (We have replaced v 2>v by — V c /2 in doing this 
transformation, which is a valid approximation since the ec- 
centricity is assumed to be small.) Obviously, if the orbit 
is circular to begin with, it stays circular. Next, in the ab- 
sence of accretion (M 2 = 0), the eccentricity will increase, 
since most of the wind material is lagging behind star 2 
{Gi,y > Gi, y > 0), as can be seen on Fig. lb of Paper I. 
Hence the second term is positive. Finally it will depend on 
the value of the accretion rate M 2 and the wind velocity at 
the point it is accreted, whether e will increase (second term 
dominates) or decrease (first term dominates). 

Applying Eq. ( |l5] ) to our case, we find that the eccen- 
tricity decreases for our set of parameters (see Fig. its 
relative variation being of the order of —dM 2 /Al 2 , i.e. a few 
percent only. Including the gravitational force due to the gas 
on the mass-losing star strengthens this conclusion. 
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4.3 Influence on spin 

The companion star may accrete spin angular momentum 
from the wind, which may alter its rotational velocity. Pack- 
ett (1981) showed that a star needs to accrete only a few 
percent of its own mass from a keplerian accretion disc to 
be spun up to the equatorial centrifugal limit, essentially be- 
cause stellar moments of inertia are generally much smaller 
than MR 2 . Less drastic effects than centrifugal breakup 
could already be detectable on a statistical basis by com- 
paring the rotational velocities of post-mass transfer stars 
with those of single stars of the same age, thus suggesting a 
further observable diagnostic of mass transfer. 

The initial spin angular momentum 5*2,0 of the star is 
defined in terms of its moment of inertia 72, o through 



&,o = £2,0 ^2,0 
h = M 2 r 2 q Rl, 



(17) 
(18) 



where ^2,0 is the initial value of I 2 . The dimensionless quan- 
tity r g is the gyration radius. The initial spin velocity of 
the star is related to its initial rotational velocity through 
^2,0 = v to t,o/R2, with R2 the stellar radius. The accretion 
of angular momentum AL will change the angular velocity 



from uj 2 ,o to tu 2 = (6*2,0 + AL)(M 2 , + AM 2 )~ 



" R 



suming that r g and R2 did not change much as a result of the 
accretion process. This value should not exceed the critical 
centrifugal angular velocity lu ci = [G(M 2 + AM 2 )/R 2 ] 1/2 , 
with G the gravitational constant. In terms of the accr eted 



spin momentum ALAM 2 1 A 1 V C 1 as defined in Sect. 3.1 
the rotation velocity v Iot after spin accretion writes 



«rot = («rot,0 + 



AL A1AM2 V.)/(1 + ^)- (19) 



AM 2 AV C R 2 rf M 2 , 



The above equation implicitly assumes that the accreted 
spin has been redistributed over the entire envelope, which 
is equivalent to assuming that the accreted matter has been 
fully mixed in the convective envelope. 



5 BARIUM STARS 

Though the simulations described in this paper could be 
applied to a variety of binary systems, they were aimed at 
gaining some insight in the origin of barium stars. In partic- 
ular they aim at checking the validity of the wind accretion 
model for the formation of barium stars developed by Boffin 

6 Jorissen (1988). The impact of our hydrodynamic results 
on the validity of that scenario will therefore be briefly dis- 
cussed in this section. 

In order to evaluate the total variation of the orbital 
elements due to the wind accretion process, Eqs. (p^)-(p^[) 
have to be integrated with Mi as the independent variable 
until the mass of the AGB primary has been reduced from 
its initial value M^o to that of its WD descendant. The total 
amount of mass lost by the primary during its evolution from 
AGB to WD, AMi, can be derived from the initial - final 
mass relationship. According to Weidemann (1984), a star 
of initial mass Mi,o will give rise to a WD of mass Mi 15 , 
such that 



MY D = 0.528 - 0.044Mi, + O.OI6M1 . 



(20) 



holds true as well for WD's in detached binary systems where 
mass loss occurred through a wind. 

To obtain changes in composition and changes in or- 
bital parameters, we start with Mi = Mi,o and integrate 
Eqs. (|lll), @, @ and @ to find M 2 , A, P and e. The 
integration is stopped when the mass-losing star reaches its 
WD mass. The total amount of mass accreted by the barium 
star, AM2, can then be used to predict the dilution factor 
T appearing in Eq. (|l^) and, hence, the pollution of the 
barium star envelope. 



5.1 Mass accretion rate 

The major result of the present paper is that the mass accre- 
tion rate obtained from our SPH simulations is significantly 
lower than the commonly used Bondi-Hoyle value. The ex- 
tent to which such a reduction affects the ability of the wind 
accretion scenario to account for the chemical peculiarities 
of barium stars is now discussed. 

As a first order approximation we assume that the ac- 
cretion efficiency remains the same during the whole AGB 
mass loss episode, i.e. /3 = 0.01 (as discussed in Sect. 3.1) 
throughout. In that case, of the ~ 2.45 Mq lost by the 3 Mq 
AGB primary [see Eq. 0.02 M Q will be accreted by 

the companion. This enriched material will be diluted in 
a ~ 1.1 Mq envelope (i.e. adopting 0.4 Mq for the core 



Iwayf 



and adopting 



mass), yielding a dilution factor T = 0.02. Assuming typi- 
cal heavy-element overabundances of the order of g = 100 in 
the AGB wind (e.g. Utsumi 1985) , we find a final overabun- 
dance log 10 / = 0.5 dex from Eq. ([[2]), quite a typical over- 
abundance factor for barium stars (see e.g. Lambert 1985). 
In the following, dilution factors T will al way s be converted 
into overabundances / by means of 
g = 100 and M 2 , C orc = 0.4 Mq. 

The accretion efficiency j3 will probably not remain con- 
stant during the wind accretion process, since it very likely 
depends upon the orbital separation (through the ratio k of 
orbital to wind velocity) and upon the masses of the two 
stars, as suggested by Eq. (|ll]) derived in the framework of 
the Bondi-Hoyle theory. It is therefore necessary to take into 
account the variations of the parameters k and fi resulting 
from th e dy namics of the mass transfer process, as outlined 
in Sect. 3.2. Unfortunately, the small number of SPH simu- 



It is assumed that this relation, derived for single WD's, 



lations performed prevents us from checking the functional 
dependence between f}, k and fj, suggested by Eq. (p"l|). Con- 
sequently, we are restricted to use Eq. ( p4| ) , yet we normalize 
it in such a way as to yield our SPH result ft = 0.01 when 
v w = 15 km s" 1 , A = 3 AU, Mi = 3 M Q and M 2 = 1.5 M Q . 

Starting with the initial parameters as listed above, the 
average (3 over the whole mass transfer process is found to 
be 0.011, i.e., slightly larger than in the first order approxi- 
mation where orbital parameters are kept constant (which is 
equivalent to fixing j3 at its initial value of 0.010). In other 
words, the accretion rate increases as the system widens 
due to the mass lost from it. This may seem rather counter- 
intuitive at first, but is easily derived from Eq. (|l(]) provided 
that w wind << v olh : P <x A~ 2 v^ h oc A~ 1/2 (M 1 + M 2 y A/2 
making use of Kepler's third law, and finally f) oc A since 
Eq. (113) predicts A cx (Mi + M 2 )^ 1 during mass transfer 
when drag effects are negligible. Note that this behaviour 
holds as long as the wind velocity is negligible with respect 
to the orbital velocity. 



© 0000 RAS, MNRAS 000, 000-000 



10 T. Theuns, H.M.J. Boffin, A. Jorissen 



For wide binaries, smaller /3 values prevail, since one can 
now neglect the orbital velocity with respect to the wind ve- 
locity to find that (3 scales as 1/A 2 according to Eq. (|f"o|). 
It is therefore expected that wind accretion will be less effi- 
cient in wider systems. The sensitivity upon k predicted by 
Eq. ( |ll| ) appears to be rather small, however, since it pre- 
dicts average f3 values of 0.011, 0.009 and 0.007 for initial 
orbital separations of 3, 4 and 5 AU, respectively, all other 
parameters being as listed above. A substantial pollution of 
the companion's envelope will still result from these values 
of j3, since the corresponding dilution factors T — 0.023, 
0.018 and 0.014, respectively, yield log 10 / = 0.52, 0.44 and 
0.38 dex. 

The threshold for producing a barium star in a binary 
system may be set to T > 0.005, yielding overabundances 
log 10 / > 0.17 dex, according to Eq. (JT^) . Below this value, 
the weakly polluted star would no longer be flagged as a 
barium star. With our new SPH scaling of Eq. (0), we 
find that the dilution factor T drops below 0.005 for sys- 
tems with final orbital periods in excess of 90 y (adopting 
the same parameters as above, except for the initial sep- 
aration). The barium star with the longest known orbital 
period, ( Cygni (P = 17.8 y, Griffin 1991), exhibits an aver- 
age overabundance of heavy elements amounting to 0.4 dex 
(Zacs 1994). Such an overabundance can easily be accounted 
for with the above prescriptions, which predict T = 0.016 
and log 10 / = 0.4 dex as required. 

We conclude that the reduction in the wind accretion 
efficiency derived from our detailed SPH simulations when 
compared to the Bondi-Hoyle prescription does not endanger 
the validity of the wind accretion scenario. In particular we 
find that wind accretion can easily account for the level of 
pollution in barium stars with orbital periods as long as 
90 y. In fact, Boffin & Zacs (1994) already suggested that a 
reduction of the accretion rate by a factor of ten with respect 
to the Bondi-Hoyle value improves the agreement of wind 
accretion predictions with observed levels of enrichement in 
barium stars. 

However, that scenario faces serious difficulties for sys- 
tems with short orbital periods, albeit on different grounds. 
Indeed, barium systems with orbital periods of a few hun- 
dred days could not stay detached when the primary star 
evolved up the AGB. Standard prescriptions for AGB evo- 
lution (e.g. Groenewegen & de Jong 1993) indicate that the 
mass-losing star will overflow its Roche lobe at some point 
during its AGB evolution in systems with P < 3000 d. Since 
the star has a deep convective envelope at that stage, the 
severe mass loss due to the Roche lobe overflow (RLOF) 
will cause its radius to expand. But the Roche lobe tends 
to shrink as a result of mass transfer from the more massive 
to the less massive component, thus leading to a runaway 
behaviour of the mass transfer process. A common envelope 
surrounding both stars is expected to result from this pro- 
cess, and the drag exerted by the envelope on the binary 
motion will cause a severe loss of orbital energy. Binary sys- 
tems with very short orbital periods (a few hours), like e.g. 
cataclysmic variables, are formed in such a way (e.g. Meyer 
& Meyer- Hofmeister 1979). Barium systems with periods of 
a few hundred days must somehow have escaped this dra- 
matic fate, though these systems are not wide enough to 
prevent an episode of RLOF. The problem of the short pe- 
riod barium systems has been illustrated by the recent anal- 



ysis of HD 121447 (Jorissen et al. 1995a), the barium star 
with the shortest orbital period (185 d). Possible solutions 
are investigated by Han et al. (1995). 



5.2 Momentum accretion 

Boffin et al. (1993) compared orbital elements of barium 
stars against those of normal G-K giants. They concluded 
that the (e,logP) diagram of barium stars could be repro- 
duced by translating to larger periods the points represent- 
ing normal giants. An increase of the orbital period is easily 
explained by mass loss from the system [see Eq. ©], but 
the change in eccentricity depends on both the accretion 
rate and the influence of the gravitational and hydrodynamic 
forces [Eq. @]. Unfortunately, the latter equation is prob- 
ably only valid for nearly circular orbits (for which we find a 
small decrease in eccentricity). Hence, we cannot draw defi- 
nite conclusions about de/dt for systems which start with a 
finite eccentricity. In addition we caution that it is very un- 
likely that the barium stars with the shortest orbital periods 
result from the widening of systems with even shorter peri- 
ods, in view of the difficulties raised by RLOF (see Sect. 5.1). 



5.3 Spin accretion 

The mild barium star HD 165141 possibly bears the signa- 
ture of spin up due to angular momentum accretion during 
the mass transfer process (see Jorissen et al. 1995b for a de- 
tailed discussion of that system). HD 165141 is a K0III Bal 
star rotating much faster (vsini — 14 km s -1 ) than normal 
K giants or barium stars. As shown by de Medeiros, Mayor 
& Simon (1992), the usini distribution of normal (single) 
K giants peaks at very low values (vsini < 1 km s _1 ) and 
decreases rapidly with increasing vsini. The median of this 
distribution amounts to 2 km s _1 , with 97% of the stars 
having vsini < 10 km s _1 . 

The rapid rotation of HD 165141 probably accounts for 
the unusually high level of chromospheric activity exhibited 
by that star as compared to other barium systems (Call 
and CIV A1550 emission, relatively hard X-rays detected 
by ROSAT). Actually, the system has typical features of a 
RS CVn system, except for the long orbital period of 5200 d. 
This long orbital period makes it very likely that wind ac- 
cretion rather than RLOF is responsible for the chemical 
peculiarities of this barium star. Our predictions for spin 
accretion presented in Sect. 4.3 may thus be applied to the 



HD 165141 system, which is too wide for tidal coupling to 
be responsible for the fast rotation rate, as is the case in 
genuine RS CVn systems. 

Adopting ALAM 2 ~ 1 A~ 1 V~ 1 ~ 0.06, as read from Ta- 
ble ^, we use i?2 = 10 R© (see below), r g — 0.3 (e.g. Tout 
& Hall 1991), Mi = 0.6 MojMa = 1.5 Mq, V c = 16 km s _1 , 
Wrot,o = 2 km s _1 in Eq. (hpl) and find that the spin up of 
HD 165141 requires a fraction AM 2 /M 2 , = 0.008 to be ac- 
creted to be consistent with its high current spin rate. The 
associated chemical pollution is consistent with that star 
displaying only mild chemical anomalies, since T = 0.011 
and log 10 / = 0.3 dex. 

The reason why HD 165141 seems to be the only bar- 
ium star exhibiting such a signature of spin accretion may 
be related to the fact that it is also the barium star with 
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the hottest (i.e. youngest) WD companion (having T c g ~ 
35000 K; Fekel et al. 1993). That situation contrasts with the 
other barium stars where the WD is usually not detected by 
the IUE satellite, indicating that it must be much cooler (i.e. 
older; see e.g. Bohm-Vitense et al. 1984). The cooling time 
scale of the WD companion of HD 165141 (twd ~ 10 7 y) is 
in fact shorter than the time spent by the present barium 
star as a giant, either on the first giant branch (assuming 
in this case that it is less massive than ~ 2M©) or in the 
He clump. It may therefore be inferred that the mass trans- 
fer must have occurred when the barium star was already a 
giant star, contrary to the situation prevailing for the ma- 
jority of barium stars that were formed as dwarfs (e.g. North 
& Duquennoy 1991; North, Berthet & Lanz 1994). For the 
dwarf stars, magnetic braking as well as the increase of mo- 
ment of inertia accompanying the radius expansion on the 
giant branch will slow the rotation down (e.g. Habets & 
Zwaan 1989), thus erasing the possible spin up resulting 
from the accretion process. Magnetic braking is probably 
still operating in HD 165141, on a time scale of the order of 
10 8 y if v Iot ~ 14 km s _1 [see Eq. (2) of Habets & Zwaan 
1989]. The slowing down of HD 165141 during the last 10 7 y, 
corresponding to the time twd elapsed since the end of the 
mass transfer process, has thus been rather moderate, in 
agreement with the large observed spin velocity. 

The current radius of R2 = 10 Kq adopted for 
HD 165141 in the above analysis corresponds to the lower 
limit derived from its 35 d photometric period and its 
vsini = 14 km s _1 rotational velocity (Fekel et al. 1993). 
With such a radius and T c g<^ 4700 — 5000 K, the star must 
currently be either in the core He- burning phase of evolution, 
or about midway on the RGB, according to the evolutionary 
tracks of Schaller et al. (1993). At the time of accretion (i.e. 
about 10 7 y ago, as given by the cooling time of the WD 
companion), the stellar radius may have been very different 
from its current value, being either much smaller (in the case 
that accretion occurred when HD 165141 was on the lower 
RGB, assuming it is still on the RGB now) or much larger (in 
the case that accretion occurred when HD 165141 was on the 
upper RGB, assuming that it reached its present core He- 
burning stage less than 10 7 y ago). It would be similar only 
in the case that HD 165141 remained in the core He-burning 
phase for the last 10 7 y (which is actually the most probable 
situation, as core He-burning lasts for at least 2.4 10 8 y in 
solar- metallicity stars of mass 2M0 or less). The fact that 
the stellar radius possibly varied in the course of stellar evo- 
lution following the accretion event, does not, however, in- 
validate the previous analysis based on Eq. (|l!]), provided 
that the spin velocity before accretion, « ro t,o, against which 
the current spin rate v is compared, is taken equal to the 
typical spin rate of normal stars in the same evolutionary 
phase, i.e. having the same current radius as HD 165141. In 
doing so, possible radius changes cancel out as they affect 
v and v rot ,o in the same way (provided that the accretion 
process did not alter the gyration radius r g ). 

The possibility that dwarf barium stars have accreted 
spin angular momentum can be investigated on the sample 
studied by North & Duquennoy (1991) and North, Berthet 
& Lanz (1994). With the exception of one star (HD 198583), 
which has a rotational velocity substantially larger than F 
dwarfs of the same temperature, the other barium dwarfs 
seem to be rotating at rates comparable to normal (non- 



barium) dwarfs. However, as in the case of barium giants, 
no WD companion could be detected for the barium dwarfs 
observed with IUE by North & Lanz (1991), implying ages 
> 2 10 s y. Since such ages are larger than the magnetic- 
braking time scale estimated above, the absence of fast ro- 
tators among dwarf barium stars is not at all unexpected. 
Ultraviolet data for HD 198583 would be of interest, in order 
to confirm our prediction of a correlation between fast rota- 
tion and the presence of a hot (i.e. young) WD companion. 



6 SUMMARY 

The SPH simulations presented in this paper model the ac- 
cretion of the wind of a mass-losing star by its companion 
in the case where the wind velocity is comparable to the 
orbital velocity. The binary motion has then a strong influ- 
ence on the accretion process, making analytic predictions 
based on the Bondi-Hoyle model suspect. In fact, we find 
that the mass accretion rate is at least ten times smaller 
than the rate predicted by the Bondi-Hoyle prescription, 
amounting to about 1% in a binary system with Mi = 3 Mq , 
M2 = 1.5 Mq, A = 3 AU and v W md = 15 km s _1 (in the case 
of an adiabatic gas with 7 = 1.5). Smaller polytropic in- 
dices lead to slightly larger accretion rates of the order of 
8%. These values probably represent upper limits, as the ac- 
cretion rate decreases with increasing numerical resolution. 

Despite the fact that the efficiency of wind accretion 
appears to have been overestimated in previous studies re- 
lying on the Bondi-Hoyle prescription, wind accretion still 
appears efficient enough to account for the chemical pecu- 
liarities exhibited by barium stars with orbital periods up 
to about 90 y. 

The accretion of transverse linear momentum, control- 
ling the eccentricity variation, appears to be negligible, at 
least in the present simulation involving a circular orbit. 

Accretion of spin angular momentum is substantial, 
with a well-developed accretion disc forming in the isother- 
mal case, contrary to a widespread belief that wind accretion 
cannot lead to the formation of such a disc. The spin up of 
the accreting star resulting from wind accretion may even be 
considered as an important diagnostic of this process. The 
mild barium star HD 165141 may be one case where such 
a spin up took place recently. Such fast-rotating post-mass 
transfer objects must be observed before the various brak- 
ing processes (among which magnetic braking) slow the star 
down again. Therefore, we predict a correlation between WD 
temperature and barium star spin rate: fast spinning barium 
stars are predicted to have hot WD companions. 
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APPENDIX 

Following the method of Huang (1956) we derive expressions 
for the change in orbital parameters due to accretion and the 
action of gravitational forces. Assuming that star 1 suffers a 
spherically symmetric mass loss, then conservation of linear 
momentum in a small interval of time dt can be written as: 

Gidi = [(Mi+dMi)(vi + dvi)-dMiVi] 



-[Mm] (21) 
G 2 dt = [(M2 + dM 2 )(v2 + dv 2 )] 

- [M 2 v 2 + 9m w v w ] , (22) 

w 

where the G's refer to the gravitational force exerted by the 
wind on the corresponding star and ^2 dm w v w denotes the 
change in momentum of the accreting star due to momentum 
accreted from the wind. The change in orbital parameters 
A, P and e follows from the change in mass, energy E and 
angular momentum h per unit reduced mass through: 

dA/A = 9(Mi + M 2 )/(Mi + M 2 ) - dE/E (23) 
dP/P = d(Mi + M 2 )/(Mi + M % ) - ^dE/E (24) 

ede/(l - e 2 ) = d(Mi + M 2 )/(Mi + Ms) - \dEjE 

-dh/h, (25) 

where E = vl /2 + vf./2 - (Mi + M 2 )/r and h = rv t , with 
r,v r (— r) and vt being respectively the radius vector, the 
radial velocity and the transverse velocity of the accreting 
star in the relative frame where the mass-losing star is at 
rest. Denoting by vi, y and v 2:V the projections of the indi- 
vidual (inertial) velocities on the y direction transverse to 
the radius vector and oriented along the orbital motion of 
star 1, one has v t = vi lV — v 2 , y - Orbit-averaged variations, 
denoted by (), can be computed from the change in the rel- 
ative transverse velocity dv t , obtained from Eqs. ( pl| ) and 
©: 

dv t — dvi sV — dv 2:V 

= -Fydt/M 2 

= Gi, y dt/Mi - G 2 , y dt/M 2 

— 'y^^dm w (v w ,y — v 2 ,y)/M 2 . (26) 

Here, as before, Gi lV and G 2>v are the projections of the 
gravitational forces on the transverse direction. With the 
above conventions, Gi tV and G 2 , y are positive, since gas is 
lagging behind the accretor with respect to its orbital mo- 
tion. Hence, the mass-losing star is accelerated by that mate- 
rial whereas the accretor is decelerated. Note that dvt refers 
to the variation of vt induced by the mass transfer (and not 
to the variation along the orbit due to keplerian motion). 

When computing the orbit averages (dE) and (dh), it 
is assumed that the mass accretion rate M 2 and the force 
F y do not depend on the phase of the binary. Clearly, this 
assumption will be poor when the eccentricity is large, so 
that the relations derived here are only valid for systems 
with small eccentricities. It should also be noted that keep- 
ing the force F y constant along the orbit is not equiva- 
lent to Huang's parametrization in terms of constant A = 
-d{M 2 v t ) /V c dM 2 = -Vt/Vc + F y dt/V c dM 2 . Indeed, A and 
F y cannot both remain constant along an eccentric orbit, 
since Vt/V c varies with orbital phase in an eccentric orbit. 
Therefore, the relations derived here are not equivalent to 
Huang's. 

Averaging over the orbital period then amounts to using 
the following mean values: 

(r) = (27) 
<r 2 ) = I4 2 (l-(l-e 2 ) 1/2 ) (28) 
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(r) = A{l + e 2 /2) (29) 
(1/r) = l/A (30) 
(1/r 2 ) = l/[A 2 (l-e 2 ) 1/2 ]. (31) 

As an example of how the calculation is done, we will com- 
pute the orbit-averaged change in kinetic energy per unit 
mass, T = ^v 2 + ^r 2 . Assuming F y (i.e. dvt) constant along 
the orbit, the orbit- averaged differential of the first term in 
T writes: 

(vtdvt) = (vt)dvt 

- (32) 
A M 2 1 ' 

where we used the average (v t ) = h(l/r) = h/A. 
Analogously, we obtain for the second term: 

(rdr) = (rd(^)) 

,.d(M 2 f) . 2 0M 2 . 
— / r — ^ -— r > 

V M 2 M 2 ' 

= -V?(l-(l-eY /2 )^ (33) 

where we assumed d(M 2 r) to be constant over an orbital 
period and used the averages in Eqs. (27-28). The sum of 
these last two results gives (dT). 

Proceeding in the same spirit, we get 

E ~ Ah+M 2 {6) 
(dh) _ (l + e 2 /2) F y dt 



h (1 - e 2 )V2 V C M 2 



(35) 



where fi denotes the potential energy per unit mass. Using 
these orbital averages in Eqs. (23-25) then leads to Eqs. 
at order e 2 . 

The change in eccentricity given by Eq. (^i|) can be un- 
derstood intuitively as follows: first recall that, for a given 
energy, the circular orbit has the maximum angular mo- 
mentum, ft follows that increasing the angular momentum 
dh > while keeping masses and energy constant, decreases 
the eccentricity, hence circularising the binary. Increasing 
the energy dE > 0, while keeping h constant, has the oppo- 
site tendency. 
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